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We study Majorana zero energy modes (MZEM) that occur in a s-wave superconducting surface, 
at the ends of a ferromagnetic (FM) chain of adatoms, in the presence of Rashba spin-orbit interac¬ 
tion (SOI) considering both non self-consistent and self-consistent superconducting order. We find 
that in the self-consistent solution the average superconducting gap function over the adatom sites 
has a discontinuous drop with increasing exchange interaction at the same critical value where the 
topological phase transition occurs. We also study the MZEM for both treatments of supercon¬ 
ducting order and hnd that the decay length is a linear function of the exchange coupling strength, 
chemical potential and superconducting order. For wider FM chains the MZEM occur at smaller 
exchange couplings and the slope of the decay length as a function of exchange coupling grows with 
chain width. Thus we suggest experimental detection of different delocalization of MZEM in chains 
of varying widths. We discuss similarities and differences between the MZEM for the two treatments 
of the superconducting order. 


I. INTRODUCTION 

Non-magnetic impurities in a conventional supercon¬ 
ductor have no pair-breaking effect and do not lead to 
qualitative change^, in contrast to the case of unconven¬ 
tional gapless superconductors where their effect is pair- 
breaking^. On the other hand magnetic impurities in 
a conventional superconductor have interesting effects^. 
A single magnetic impurity coupled to the conduction 
electron spin density gives rise to a local bound state, 
known as Yu-Shiba-Rusinov (YSR) state, that is created 
by breaking a Cooper pair and capturing an electron with 
the appropriate spin component. As the exchange 
coupling J increases, a first order quantum phase tran¬ 
sition (QPT) occurs, at which point the gap function 
has a TT shift and the magnetization of the conduction 
electrons jumps from zero to 1/2. This phase transition 
has been detected experimentally^^”^^ and is signalled in 
various quantities including quantum information signa- 
tures^^’^^. As it turns out^“^ the YSR states come in 
pairs that lay inside the gap (one at positive and one 
at negative energy) and tend to lower energies as the 
coupling grows. At the quantum critical point there is 
a level crossing, such that the bound state has a small 
but finite energy. Increasing the number of impurities, 
more states appear inside the gap (two per impurity) 
and as the coupling increases, a series of quantum phase 
transitions occurs and the magnetization of the conduc¬ 
tion electrons changes in increasing plateaus. Interesting 
interference phenomena occur as the number of impuri¬ 
ties grows and extended states appear along the chain of 
impurities placed on the two-dimensional superconduc- 
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A quantum treatment of impurities with spin 1/2 de¬ 
scribes the competition between the screening of the im¬ 


purity spin and the superconducting correlations, which 
lead to a similar type of QPT between the spin-doublet 
and spin-singlet state for the unscreened and screened im¬ 
purity spin, respectively^^“^^. Real magnetic impurities 
however, have higher spin, where magnetic anisotropy 
plays an important role^^. The effects of magnetic 
anysotropy on the subgap excitations induced by quan¬ 
tum impurities in conventional superconductor^^ were re¬ 
cently measured^^ in Manganese pht halo cyanine (MnPc) 
on top of superconductor Pb, where Mn has high spin 
3/2. Different possible ways of adsorption of MnPc 
molecules enabled the precise spectroscopy of YSR states 
across the QPT using scanning tunnelling microscope 
(STM). 

Increasing the number of magnetic impurities eventu¬ 
ally leads to destruction of superconductivity. If the mag¬ 
netic impurities have arbitrary orientations and locations 
on the superconductor the pair breaking effect leads to 
gapless superconductivity and eventually destruction of 
superconductivity^"^ for small concentration of impurities 
of the order of a few percent. A higher robustness of 
superconductivity to the increasing number of magnetic 
impurities has been found if they are correlated, partic¬ 
ularly if their locations are not random but organized in 
some patterns^^. These regular distributions allow quite 
high concentrations of impurities without destruction of 
superconductivity. 

In 2014 an experiment was reported^^ in a system of 
a chain of magnetic adatoms (Fe) placed on top of a 
two-dimensional conventional superconductor (Pb). Lo¬ 
calized zero energy modes were detected at the edges 
of adatom chain using STM, which seems to provide 
evidence of MZEM. These zero energy modes have at¬ 
tracted considerable attention recently^^“^^ because of 
their non-abelian statistics, which can be used for quan- 
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turn information storage and manipulation^^^^^. A key 
model exhibiting such topological properties is the one¬ 
dimensional Kitaev model of spinless fermions with p- 
wave superconducting pairing^^. More feasible models 
based on a conventional s-wave superconductor in prox¬ 
imity to either topological insulator^^ or semiconducting 
nanowires in the presence of a Zeeman field^^’^^ greatly 
enhanced the interest in the field. The later proposal ap¬ 
pears to be experimentaly realized in a zero-bias conduc¬ 
tance peak measurements in a semiconducting nanowire 
on top of a conventional superconductor^^’"^^. Although 
all these experiments support the existence of MZEM 
they are for now inconclusive^^”^^. 

In the experiments^, the zero energy states are sharply 
located at the edges. Complementing the experiment 
was a theoretical study of the system, which in its sim¬ 
plest form considers 2D superconducting surface with 
nearest neighbour hopping, chemical potential, exchange 
interaction with the ferromagnetically aligned magnetic 
adatoms, Rashba SOI and proximity induced s-wave su¬ 
perconductivity. As parameters of the model change, the 
system evolves from trivial to topological phase, as deter¬ 
mined by the non self-consistent solution of the problem, 
justifiable by a proximity induced superconducting order 
or by neglecting the change of the gap function and con¬ 
sequently of other physical quantities in the vicinity of 
the magnetic adatoms. 

Moreover, there are many theoretical studies show¬ 
ing that regularly positioned chains of adatoms give 
rise to MZEM. These were given for the cases of ran- 
dom"^^, spiral^^“^^, AEM^^’^^ and order¬ 

ings of the magnetic adatoms. A stable spiral mag¬ 
netic structure of the ID adatom chain is established 
via the Ruderman-Kittel-Kasuya-Yosida (RKKY) inter- 
action^^“^^ mediated by the electrons in the underlay¬ 
ing superconductor, which has to be treated fully self- 
consistently^^'^^. A more realistic 2D model of a super¬ 
conducting surface coupled to a ID chain of magnetic 
moments shows^^’^^, that besides spiral ordering, there 
are regions in the phase diagram where the magnetic mo¬ 
ments arrange either EM or AEM. These however, can 
in the presence of SOI, still host MZEM which occur 
in multiple topological phases^^’^^ even in the presence 
of periodic modulation of superconducting order param¬ 
eter^^. Recently magnetic chains and the competition 
between MZEM with the zero energy modes at the bor¬ 
der of a triplet superconductor were also studied^^. The 
importance and effects of the self-consistent solution on 
the order parameter in triplet superconductor were also 
discussed very recently^^. 

In this work we are interested in comparing the results 
obtained from the non self-consistent solution with those 
obtained by treating the superconducting gap function 
self-consistently. The case of a few impurities has shown 
the relevance of a full treatmentand we aim to see 
how the results change as we pursue a more detailed solu¬ 
tion. As was already shown previously for spiral magnetic 
chains^^’^^ qualitatevely both treatments give similar re¬ 


sults and here we confirm such results for the case of 
EM chain in the presence of Rashba SOI, with notable 
quantitative shift of the regions with MZEM in chemical 
potential. We study also Majorana decay lengths and 
find that as a function of exchange interaction, chemical 
potential and the superconducting gap function (effective 
attractive interaction for the self-consistent superconduc¬ 
tivity) the decay length is a linear function with positive, 
zero and negative slope, respectively. By considering EM 
chains of various widths, we show that the wider chains 
correspond to a narrow chain with increased effective ex¬ 
change interaction. Due to the linear dependence of the 
Majorana decay length on exchange interaction differ¬ 
ent delocalizations of MZEM in the STM measurements 
should be detectable on chains of various widths. 

An important difference between the two treatments of 
superconductivity is that in the self-consistent not only 
topological phase transition but also a series of quan¬ 
tum phase transitions occur. The later are due to vari¬ 
ous ingap level crossings and can be observed in changes 
of total magnetization of the electron spin density, local 
spin density, local density of states, various quantum in¬ 
formation measures and the local order parameter. In 
this work we focus on the later and show that the local 
superconducting gap on the adatom sites experiences a 
discontinuous drop at the exchange coupling where the 
topological phase transition occurs. 

The paper is organized as follows: in section II we in¬ 
troduce the considered model. In section III we present 
solutions of the model for both treatments of the super¬ 
conducting order parameter. We study energy levels in 
subsection III A, discuss topological invariant in subsec¬ 
tion IIIB study the MZEM wavefunctions in subsection 
III C and consider the case of wider chains and the Ma¬ 
jorana decay lengths in subsection HID. In section IV 
we focus on the self-consistent treatment of supercon¬ 
ductivity and study local order parameter. We present 
our conclusions in section V. 


II. MODEL 

We consider a 2D s-wave superconducting surface with 
adsorbed magnetic adatoms, which we describe as classi¬ 
cal spins and assume their EM alignment. The classical 
spins are placed along the line of length L and width w 
and oriented along the normal to the surface. The SC 
surface is a 2D square lattice as shown in Eig.l and the 
Hamiltonian of the system considered is 

H = —t ^is^js ~ (1) 

{i,j),s i,s 

+ Y {4+Sx,s<^L>Ci,s' - c4sy,s^Bs'Ci,B' + h.C.) 
i,s,s' 

+ + A^c^cit) - Y Ji4,B^Bs'Ci,B', 

i i,s,s' 
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where the first term describes the electron hopping be¬ 
tween two adjacent lattice sites, the second term is the 
chemical potential /i, the third is the Rashba SOI, forth 
term is the superconducting s-wave pairing, with site 
dependent order parameter and the last term is ex¬ 
change interaction with the magnetic adatoms located at 
specific sites. We consider = J at these adatom sites. 
Sx{Sy) are vectors to the next neighbouring sides along x 
(y) direction, are Pauli spin matrices with a = x^y^z 
and 5 , 5 ' =t,Note that the indices = 1,...,A , 
where N is the total number of sites. 
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FIG. 1. Systems considered: Example of a 2D square lattice 
of size Nx X Ny, with the magnetic adatom chain (black) of 
length L and width le = 1. 
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where f{i) = f{i + 1) i f{i — 1), for I = x, y, i.e. dis- 
placement along x and y and coefficients f{i) are the par¬ 
ticle and hole coefficients Un{i, s)^Vn{i, s). Unless stated 
otherwise we set the hopping to unity, i.e. t = 1 and 
express all the energies in this unit. 

We solve BdG equations by diagonalizing a (4A) x 
(4A) matrix, with the gap either constant, which we re¬ 
fer to as non self-consistent solution or calculated self- 
consistently, imposing at each iteration that 

\ = V{{ci^Cii) - {ci^Ci^))/2, (5) 

where V is an effective attractive interaction between the 
electrons. The superconducting gap function is calcu¬ 
lated via 

Ai = r ^ {1/2-fn) [un{i, t)] , 

n,Sn<OOD 

(6) 

where cji:) is Debye frequency, the sum runs over all the 
positive energy eigenstates with values smaller than 
fn = + 1] is the Fermi function and T is the 

temperature. In this work we set ujd = 2, T = 0.001 
and find that calculated A^ has in all cases considered 
negligible imaginary part. The dimension of the systems 
considered are Nx > L where the chain length L is up to 


We diagonalize the Hamiltonian using the Bogoliubov 
transformation 


n 


where n is the complete set of states, Un and Vn are re¬ 
lated to the eigenfunctions of the Hamiltonian (1) and 
dUdn) are the new fermionic quasiparticle creation (anni¬ 
hilation) operators. We calculate the coefficients Un{i^s) 
and Vn{i^s) by solving Bogoliubov-de Gennes (BdG) 
equations^^, which can be written as 

(3) 


with the 4A dimensional vector 




/Un{iA)\ 
V'^n (b T) / 


and the matrix 1-L at site i is 


+ ir]y ) 

0 

A* 


0 


\ 


Knt + vt) + + Jij 


(4) 


180 sites and the lateral dimension Ny is between 11 and 
31. 


III. MAJORANA ZERO ENERGY MODES 
A. Energy levels 

First let us focus on the single-particle energy levels of 
eq. (1) for self-consistent treatment of superconductiv¬ 
ity. In the absence of magnetic adatoms (J = 0) and SOI 
(a = 0), the energy states in a clean homogeneous s-wave 
superconductor are gapped, with the gap value of A. By 
adding a single magnetic impurity, two YSR local bound 
states appear within the gap^“^, one at positive and one 
at a symmetric negative energy. If the exchange cou¬ 
pling J, between the impurity spin and the conduction 
electron spin density is weak enough, the problem may 
be treated classically^^ and the impurity spin acts as a lo¬ 
cal magnetic field that interacts with the superconductor 
through the Zeeman effect. In this work we consider such 
case. Note also that J is actually a product of the ex¬ 
change coupling and the magnitude of the impurity spin 
S. As the exchange coupling J increases, the values of 
the two levels drop towards zero until a critical value of 










4 


J is reached, where a QPT occurs. With the further in¬ 
crease of J the levels start to repel and eventually merge 
into the gapped states. 

In the case of Ni impurities, there are 2Ni YSR states 
inside the gap for small J. Similar to the case of a single 
impurity, the energy values of these states approach zero 
by increasing the exchange coupling and in the case of 
an ordered array they form minibands^^. By increasing 
the exchange coupling several level crossings occur (all 
at finite energy) accompanied by a series of first order 
QPTs. 

In Fig. 2 are shown the lowest + l energies as a func¬ 
tion of exchange coupling for the chain of length L = Ni 
and width re = 1, where the local spins are FM aligned. 
Increasing the exchange coupling from low value a YSR 
miniband approaches zero energy and at some critical 
exchange coupling the lowest state becomes practically 
zero. By increasing the exchange coupling the lowest en¬ 
ergy state remains zero, while the energy of the first ex¬ 
cited state increases to a value of about 1/10 of the super¬ 
conducting gap. By further increasing J the first excited 
state, together with the whole miniband returns to the 
zero energy state and the minigap closes. Still increasing 
the exchange coupling the whole miniband moves toward 
the superconducting gap, which has all along remained 
nearly constant, as will be discussed in more detail be¬ 
low. Also shown in Fig. 2 is similar behaviour of the 
lowest 1/ + 1 energies at fixed exchange coupling as a 
function of chemical potential. There are several mini¬ 
gap closings and reopenings and two regions with the 
lowest energy state remaining in vicinity of zero, while 
the minigap increases. At fixed exchange interaction J 
the superconducting gap reduces with increasing chemi¬ 
cal potential. We also note that the superconducting gap 
function is in the case of self-consistent superconductivity 
site-dependent and changes its average value by chang¬ 
ing either chemical potential, SOI or effective attractive 
interaction. Similar effects were already observed in var¬ 
ious magnetic and superconducting hybrid systems^^”^^ 
and studied theoretically in magnetic chains on top of a 
superconducting surface^^’^^ in the absence of SOI. 

The results for non self-consistent treatment of super¬ 


conductivity are qualitatively similar for lowest energy 
states, but quite different when considering lowest L + 1 
energy levels. The lowest L energy states are within the 
superconducting gap A only in vicinity of zero chemical 
potential. Away from zero chemical potential no clear 
superconducting gap is observed. 

Next we focus on the lowest single-particle state en¬ 
ergy. Examples of contour plots as a function of chemi¬ 
cal potential /i and exchange coupling J for fixed values 
of superconducting gap A (effective attractive interac¬ 
tion V) and SOI a are presented in Fig. 3 for non self- 
consistent (self-consistent) superconductivity. Note that 
the energy scale used is logarithmic, with the cutoff at 
the energy 10“^. Blue (red) areas correspond to lower 
(higher) energies. The low energy area is mainly concen¬ 
trated in the vicinity of the chemical potential /i ^ 4, 
/i ^ 1 for non self-consistent, self-consistent supercon¬ 
ductivity, respectively. In both cases the low energy area 
appears at some minimal exchange coupling, which in¬ 
creases with increasing superconducting gap (effective at¬ 
tractive interaction V) in the case of non self-consistent 
(self-consistent) treatment of superconductivity. 


B. Topological invariant 


To characterize the zero energy states we calculate the 
Majorana number^^ M, which shows the presence (M = 
— 1) or absence (M = 1) of MZEM at the edges of chains 
in the system. Eor translationally invariant system the 
Majorana number is given as 


M = Sign[pi{Ap)Pf{Aa) 


(7) 


where Pf(A) denotes the Pfafhan of a matrix A and 
Ap (Aa) is the Hamiltonian with periodic (anti-periodic) 
boundary conditions, rewritten in Majorana basis. Intro¬ 
ducing Majorana fermion operators a 2 i-i,s = 

Ci 2 i,s = — 4 s)’ which holds aj ^ = ai^s and 

cii,sCij,s' cij,s'Cii,s = the Hamiltonian (1) is 


^ <^ 2 i-l,s<^ 2 i,s + t [^i2i-l,s<^2(i+(5),s + <^2(i+(5)-l,s<^2i,s] ~ 

i s s,6=6x,6y 

— C^[a2(i-^5x)-l,t^2i,i + 0.2i-l,iCi2(i-^Sx),t “ <^2(i+(5cc)-14<^2i,t “ <^22-l,t<^2(i+(5cc)q “ 

- (^2{i+6y)-l,t^2i-l,i - <^2(i+(5y)-14<^2i-l,t “ (^2{i+6y),t^2i,i “ a2(i+<5y)4<^2i,t] + 

+ Re(A^) [a 2 i-iqa 2 iq — Ci2i-i,iCL2i,t] + Iia(A^) [a 2 i,^a 2 iq — CL 2 i-i,^Ci 2 i-i,i] + 

+ Ji \ci2i-l,iCl2i,i — ^i2i-l,t^2i,t 



( 8 ) 


Let us first consider a ID chain {Ny = 1), which 
is just the discrete case of a single-channel quantum 


wire38’39 Taking A as homogeneous and using the 
Eourier transform a 2 j-i,s = 1/ 
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FIG. 2. Energies of the lowest 2(L + 1) single-particle energy states as a function of exchange interaction J and the chemical 
potential ji for the self-consistent superconductivity. Lowest energy level is shown with full red line. The second and L-th 
single-particle energy levels are given with full black line and define the minigap in the areas where the lowest level has zero 
energy. The YSR states, except the lowest one, are shown in the shaded area, which defines the ingap miniband. The L + l-th 
state is given by blue line and when well separated from the L-th state, shows the behaviour of the superconducting gap. The 
parameters used are L = 90, Ny — 11, ji — 0.6, J = 2, a = 0.3 and V — A. 



FIG. 3. Gontour diagrams of the lowest single-particle state energy as functions of chemical potential fi and exchange coupling 
J for a non-self consistent and self-consistent treatments are shown in left and right panel, respectively. Note the logarithmic 
scale used, with the cutoff at exponent —6. The parameters used are L = 90, Ny = 11 and a = 0.3. 


and a 2 j^s = Exp[—i/cj]a/c^ 2 ,s the Hamil¬ 
tonian (8) can be rewritten in the basis = 

(9) 

k 

where a skew-symmetric matrix A{k) has the following 
nonzero elements 

= Sk^ J, = A + 2ia sm{k), 

A 2 , 3 {k) = -A - 2iasm{k), A 2 , 4 (^) = Sk - J, (10) 

with Sk = —— 2tcos{k). The Majorana number can 
be calculated using eq. (7), where Ap = H(0) and Aa = 
H(7r), thus we find M = — 1 for 

yA2 + {fi- 2tf <j< yA2 + {ij, + 2ty. (ii) 


Note that the topologically non-trivial region is for ID 
case independent on the SOI. As is well known"^^, the 
ID chain in constant magnetic field (here J) and SOI is 
equivalent to a ID chain in spiral magnetic field, which 
is evident also by comparing the matrix A{k) to the one 
given in Ref."^^. 

The equivalence no longer holds for the 2D system, 
since the SOI considered in Hamiltonian (1) is homoge¬ 
neous throughout the system, whereas the spiral mag¬ 
netic ordering is limited only to adatom chain. For Ny > 
1 the topologically non-trivial region becomes also SOI 
dependent. In Fig. 4 we show the lowest single-particle 
energy states as functions of the considered model pa¬ 
rameters. Note the transitions of the lowest energy states 
from finite to zero values when the gaps close and reopen. 
When L = N^, the occurence of zero energy states co¬ 
incides with the negative Majorana number, which we 
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FIG. 4. The 40 lowest energy states as a function of model parameters: chemical potential // , exchange coupling J, order 
parameter A and spin-orbit interaction a. Unless given in the panels, the parameters used are L = Nx = 90, Ny = 11, 
= 5, J = 2,0; = 0.3 and A = 0.3. Shaded area denotes topological regions where M = —1 and the localized MZEM occur 
at the edges of adatom chain. A slight disagreement at the gap closings as a function of A is a hnite size effect, which we 
conhrmed by considering longer chain. 


calculate using the algorithm developed in Ref.^^. Note 
that in some cases the Majorana number is negative, but 
the minigap is closed, thus in these cases the lowest en¬ 
ergy states are not MZEM. 


C. Wave functions 

To study the zero energy modes in more detail we an¬ 
alyze their wave-functions. The wave-function is 

given as 

|^(i)|^ = ^K*,5)|2 + Kz,s)|2 (12) 

s 

where v are the particle and hole coefficients, respec¬ 
tively and the sum runs over spin s =t 5 i- Majorana 
bound states appear at the end of the EM chain, as can 
be seen in Fig. 5, middle panels. In case of Majorana 
states the particle and hole coefficients are related due 
to the particle-hole symmetry^^ and the two end state 
peaks are separated by taking the following linear com¬ 
binations, 

+ |u(i, J,) ^J,) I , (13) 

as shown in left and right panels of Fig. 5. This way 
the separate end states are obtained directly from the 


particle and hole coefficients calculated from the BdG 
Hamiltonian (4). 



FIG. 5. Majorana zero energy states: middle panels show 
|T|^ and left (right) panels linear combinations of particle 
and hole coefficients 7 l {'Jr) as given by eq. (13) for the cases 
Nx = 200, L = 180, Ny = 15, o = A = 0.1, ji — A and 
J = 0.7 (upper panels) and J = 1.1 (lower panels). 

Typical evolution of MZEM as a function of exchange 
coupling is presented in Fig. 6, where we show the wave- 
function |ff^(i)P and limit to i only at the sites of the 
adatom chain. Note that before the closing of the gap 
(small exchange coupling) the wave-function is spread all 
along the chain, corresponding to the YSR state, whereas 
at the transition to topological phase, the wave-functions 
become strongly localized at the ends of the chain. Unlike 
in ID systems^^’^^’^^ the transition here is more gradual. 
By further increasing the exchange coupling, the decay of 
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the wave-function becomes larger and the energy of the 
lowest energy state increases (inset in Fig. 6). 




FIG. 6. Evolution of the MZEM from YSR state: both panels 
show 1^1^ on the adatom sites and different states are dis¬ 
placed for better presentation, indicated also by thin dashed 
lines. YSR state is the lowest, black curve in lower panel, 
for which the exchange coupling J — 0.605. Second curve 
from bottom has J — 0.635 and between the following neigh¬ 
bouring curves there is difference in the exchange coupling of 
0.01. J grows from bottom to top, so that the top curve with 
localized end state has J — 0.685. Upper panel shows delo¬ 
calization of end states for increased exchange couplings. The 
y-axis scale is the same in upper and lower panel and given by 
the arrow. The logarithms of the two lowest energy states as 
a function of exchange coupling are shown in the inset. Note 
that the ground state energy (thick) drops close to zero and - 
due to enhanced overlap of the end states - grows with further 
increase of the exchange coupling. Other parameters used are 
= 200, L = 180, = 15, a = A = 0.1, /x = 4. 


D. Wider chains 


ing it analytically with mean-field theoryAs in the 
theoretical treatment supplementing experiment and ex¬ 
periment itself^^ we find here strongly localized MZEM. 
Due to the evolution of highly localized state from the 
YSR state, the decay length first decreases and we note 
that in this transition region the fitting function poorly 
describes the actual envelope function. After the tran¬ 
sition, the decay length is a linear function of exchange 
coupling (as can be seen in Fig. 7) and roughly inde¬ 
pendent on the chemical potential. With the exception 
of the transition regions, the decay length shows linear 
behaviour with negative slope also as a function of the 
superconducting gap function (effective attractive inter¬ 
action) for the non self-consistent (self-consistent) treat¬ 
ment of superconductivity. 

We consider also cases with various chain widths w. 
Note that the length dependence was studied previ¬ 
ously^^ and also various widths were considered but here 
we systematically study the width effect on Majorana 
decay lengths and show its implications for experiment. 
The general observation in this case is that wider chains 
correspond to a narrow chain {w = 1) with increased 
effective exchange coupling. This can be seen in the Ma¬ 
jorana decay lengths in Fig. 7. The wider the chain, 
the MZEM appear at smaller exchange coupling, which 
is also seen in the plot of the exchange coupling where 
the minimum decay length is shown for different chain 
widths. This minimum is the same for all the widths 
considered and has for the chosen parameters value ^ 4. 
Besides the shift of the decay length minimum also the 
slope of the decay length as the function of the exchange 
coupling increases with the increasing chain width. 



EIG. 7. Majorana decay lengths as a function of exchange 
coupling for various chain widths w are shown in left panel. 
Right panels show the exchange coupling where the mini¬ 
mum of the decay length is reached (upper panel) for dif¬ 
ferent widths. Lower right panel shows the growth of k with 
increasing chain widths, where k is the slope of the linear fit 
of ^(J) = /cJ + ^o- System dimensions are Nx — 200, L = 180, 
Ny — 15 {Ny — 16 for w — 2) and o = A = 0.1. 


We study also the wave-function decay length which 
has been modeled^^’^^ as |ff^(i)P ^ Strong 

localization of MZEM observed in experiments^ was ex¬ 
plained recently by modeling system as a linear chain of 
Anderson impurities on top of superconductor and solv- 


From these results we see that at fixed exchange cou¬ 
pling the wider chains have longer decay lengths (in case 
they are in the topological regime with the MZEM). Thus 
we suggest that the delocalization of the edge modes for 
wider adatom chains could be observed in the experi¬ 
ments. 



































IV. QUANTUM PHASE TRANSITIONS 

First let us review known results for few magnetic im¬ 
purities. In the case of a single magnetic impurity treated 
classically there is a first order QPT as a function of the 
exchange coupling, which was first pointed out by Saku- 
rai^. The transition corresponds to a level crossing where 
the superconductor becomes thermodynamically unsta¬ 
ble against spontaneous creation of localized quasiparti¬ 
cle excitation (with spin opposite to the local magnetic 
field due to impurity), which occurs for both non-self 
consistent and self-consistent (mean-field) treatment of 
superconductivity^. In the latter at the QPT both local 
superconducting gap function and the local spin density 
have discontinous jumps as can be seen in Fig. 3 of Ref.^. 
A jump in the local superconducting order is seen also 
in the full quantum treatment^^. In that case the QPT 
occurs only for antiferromagnetic coupling (J > 0)^^ at 
Tx/A 0.3, where Tk and A are the Kondo temper¬ 
ature and superconducting gap, respectively. Recently 
the spectral properties of YSR states at finite tempera¬ 
tures were studied using the numerical renormalization 
group . 

The case of two magnetic (classicall) impurities ex¬ 
hibits quantum interference effects, which depend on the 
interimpurity distance and the relative angle between 
the impurities^^. In contrast to the single impurity case 
the two impurity system can be driven through multi¬ 
ple QPTs, which may or may not be accompanied by a 
TT phase shift of the onsite superconducting gap^^. At 
QPTs the local superconducting gap function changes 
discontinuously and the spin polarization changes for 
±1/2. The gap function renormalization in the system 
of two impurities was recently revisited^^. Similarly to a 
pair of magnetic impurities also ordered magnetic chains 
undergo a series of first order QPTs as the exchange in¬ 
teraction of magnetic moments is increased^’^^’^^. These 
phase transitions are induced by numerous level crossings 
between the ingap states and can be observed in changes 
of total magnetization of the electron spin density, local 
spin density, local density of states, various quantum in¬ 
formation measures and the local order parameter. In 
this work we focus on the latter. 



FIG. 8. Typical superconducting gap function after conver¬ 
gence has a finite positive value away from adatom sites and 
a drop below the adatom chain, where the value of can 
be either positive or negative. In the latter case system has 
undergone a QPT. For the case shown Nx > L. 


In this section we consider only self-consistent treat¬ 
ment of superconductivity. The typical superconducting 
gap function profile after convergence is shown in Fig. 8. 
It is usually reduced to a small absolute value, which can 
be either positive or negative on the magnetic adatom 
sites and assumes a finite positive value away from them. 
Similar superconducting gap profile was reported also in 
the study of spiral magnetic chain on top of an s-wave su¬ 
perconducting surface^^, although in that work the QPTs 
were not discussed. 

To give a more detailed analysis we calculate two aver¬ 
age values of the gap function, first the average over all 
sites in the system. A, given as 

A = l/Ar^A,, (14) 

i 

and second the average over the adatom sites, AacI, cal¬ 
culated as 

AAd = Y. 

iCAd 

The overall average of the gap function depends very 
little on the exchange coupling when the chemical po¬ 
tential is kept fixed and its maximum value at /i = 
0 is roughly proportional to the attractive interaction 
strength V. This value decreases with the increasing 
chemical potential until it vanishes and the superconduc¬ 
tivity is destroyed at the chemical potential ju ^ 4. This 
behaviour of A is seen also by the L ± 1-st energy state 
in Fig. 2 in the parameter ranges where the ingap states 
are away from the superconducting gap. As a function 
of spin-orbit interaction the gap has the largest values 
at small chemical potentials, which decrease towards 0 
value with increasing a. For larger chemical potential 
the value of the gap at small a diminishes to 0 at /i ^ 4 
and grows to some finite value after some critical a which 
increases with increasing fi. 

Although A itself is very weakly dependent on the 
exchange coupling, AAd shows strikingly different be¬ 
haviour, as can be seen in Fig. 9. At small exchange 
coupling AAd A and it slowly decreases with increasing 
J. At some exchange coupling the local gap function 
suddenly drops. This occurs at ^ 1.6 for all differrent 
chemical potentials and attractive electron interactions 
considered. The slope change is greater at smaller chem¬ 
ical potentials and it occurs at smaller for smaller 
chemical potentials and smaller attractive interactions. 
At higher chemical potentials (/i ^ 1.8 for U = 4 and 
/i ^ 3 for y = 5.6) AAd changes smoothly and there is 
no abrupt slope change. As exchange coupling is further 
increased AAd passes zero value and becomes negative 
(tt shift). This critical exchange coupling Jc, where the 
supercondicting gap at the impurity sites changes sign is 
between 1.6 < Jc < 2 for wide range of chemical poten¬ 
tials and attractive electron interactions. By increasing 
the attractive electron interaction the points and Jc 
move closer. The tt shift is also observed by changing the 
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FIG. 9. Average superconducting gap function at adatom sites for various attractive interactions. The QPTs start in the region 
1.3 < J < 1.7 for /i < 4, where the slope of the gap function changes abruptly. For the exchange coupling line cuts the chemical 
potential for the highest curve (blue) is /x = 0 and the following curves are in chemical potential steps 0.3. Other parameters 
used are L = 90, = 11 and a = 0.3. 


SOI. The critical SOI where the transition occurs depen- 
dends on the chemical potential and is for /x ^ 0 at some 
finite value, which decreases to 0 at /x ^ 4. We observe 
no TT shift as a function of attractive interaction V. 



FIG. 10. Average superconducting gap function at adatom 
sites as a function of exchange coupling for various chain 
lengths L. A Ad has a drop where Majorana number changes 
sign (not shown). Lines are guide to the eyes. The parameters 
used are Ny = 11, ju = 1, a = 0.3 and V = 4. 

Closer inspection of the region in vicinity of Jd reveals 
that the superconducting gap function on the adatom 
sites has a finite drop as shown in Fig. 10. This drop is 
accompanied with the single particle gap closing and the 
occurence of MZEM. Consistently the Majorana number 
calculated using (7) becomes negative. Remarkably al¬ 
ready the case of L = 10 gives both a change in Majorana 
number and a drop in the superconducting gap, which 
happens practically at the same value of exchange cou¬ 
pling as for longer chains. However, the lowest single par¬ 
ticle energy is in this case ^ 0.03. Note that even shorter 
chain of length L = 5 does not show a drop in the lo¬ 
cal superconducting gap function, however the Majorana 
number also for this case changes sign (at J = 1.628). 
The size of the drop of superconducting gap is a bit larger 
for smaller system sizes and is the same for longer chains. 

The coincidence of the drop in local superconducting 


gap function and the change of sign in Majorana num¬ 
ber is perfect in the whole range of chemical potentials 
and for various attractive interactions considered where 
MZEM occur. As can be seen from Eig. 3 this happens 
around chemical potential /x ^ 1. Thus around /x ^ 0 
there are no MZEM, but as soon as they appear we find 
Jtpt = Jd- By increasing chemical potential the abrupt 
decrease of adatom gap function becomes smooth and 
the chemical potential where this occurs coincides with 
the disapperance of MZEM in the phase diagram. Since 
Jc> Jdi the superconducting gap function on the adatom 
sites changes sign already within the area of MZEM, thus 
the latter exist both at A Ad > 0 and A Ad < 0- 

V. DISCUSSION AND CONCLUSIONS 

The self-consistent treatment of superconductivity was 
already addressed in systems of magnetic chains on 
top of conventional superconductors both in the con¬ 
text of QPTs^^’^^ and MZEM"^^’^^. In the latter it was 
shown that qualitatively both non self-consistent and 
self-consistent treatment of superconductivity give very 
similar results, that is in both treatments MZEM are 
found in the phase diagrams. Here we confirm such re¬ 
sults for a 2D s-wave superconducting surface with a EM 
chain of adatoms, in the presence of Rashba SOI. We have 
shown that quantitatively the areas where MZEM occur 
shift in the phase diagrams, most notable is the shift in 
the chemical potential (seen in Eig. 3), where in the case 
of non self-consistent (self-consistent) superconductivity 
the MZEM occur at higher (lower) chemical potential, 
i.e. /X ^ 4 (/X ^ 1). In the self-consistent treatment the 
shift is due to the decrease of the superconducting gap 
with increasing chemical potential (seen in Eig. 2). Note 
however that whereas the bulk superconducting gap A 
is relatively high, the local superconducting gap on the 
adatom sites AAd is small (and even changes sign) in the 
MZEM region. 

By a systematic study of Majorana decay lengths in 
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wide chains (as presented in Fig. 7) we have shown that 
wider chains correspond to a narrow chain with increased 
effective exchange coupling. Due to the linear depen¬ 
dence of the Majorana decay lengths on the exchange 
coupling that we found, the wider chains have larger 
decay lengths (at fixed exchange coupling). From this 
analysis we suggest that STM experiments performed on 
wider chains should observe greater delocalization of the 
zero energy modes. 

To study the QPTs we have focused on the local su¬ 
perconducting gap function AacI (Fig- 9 and 10). As a 
function of exchange coupling J it decreases gradually 
with increasing J until a certain value Jd, where there is 
a drop, which is already present at short adatom chain 
L ^ 10 and is only weakly dependent on its length. By 
further increasing the exchange interaction the supercon¬ 


ducting gap AacI changes sign (tt shift). In this work 
we have shown that the value Jd coincides with the oc¬ 
curence of MZEM at the ends of the chain, thus there is a 
connection between the topological phase transition and 
the quantum phase transition. It is interesting to note 
that there is an analogy between the coincidence of a tt 
drop in superconducting gap at the site of single impu¬ 
rity and the QPT and the coincidence of the topological 
phase transition and the drop of local superconducting 
gap function at the adatom chain. 
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